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ABSTRACT 

Thirty eight long-period (P > 10 days), apparently contact binary stars discovered 
by the OGLE-II project in the SMC appear to be Lyrae-type systems with ellipsoidal 
variations of the cool components dominating over eclipse effects in the systemic light 
variations and in the total luminosity. A new period-luminosity-color (PLC) relation 
has been established for these systems; we call it the PLC-/? relation, to distinguish it 
from the Cepheid relation. Two versions of the PLC-/? relation - based on the (B — V)q 
or (V — I)o color indices - have been calibrated for 33 systems with (V — I)o > 0.25 
spanning the orbital period range of 11 to 181 days. The relations can provide 
maximum- light, absolute- magnitude estimates accurate to eMy ~ 0.35 mag. within 
the approximate range —3 < My < +1. In terms of their number in the SMC, the 
long-period /? Lyrae-type binaries are about 50 times less common than the Cepheids. 
Nevertheless, their large luminosities coupled with continuous light variations make 
these binaries very easy to spot in nearby galaxies, so that the PLC-/? relation can 
offer an auxiliary and entirely independent method of distance determination to nearby 
stellar systems rich in massive stars. The sample of the long-period /? Lyrae systems 
in the SMC analyzed in this paper is currently the best defined and uniform known 
sequence of such binaries. 



1. INTRODUCTION 



This is the fifth paper based on the OGLE data, but the first one based on the data from 
the second phase of the OGLE project (OGLE-II, Udalski et al. 1997| ) utilizing the new 1.3 meter 
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Polish telescope on Las Campanas. The data are freely available over the computer networks^]. 

Analysis of large numbers of eclipsing binary stars discovered by the OGLE-I project in the 
direction of Baade's Window (BW) revealed that among almost one thousand binaries, about two 
thirds were contact systems of the W UMa-type (Rucinski 1997a, Rucinski 1997b = Papers I and 



II). The direction of the BW is towards the Galactic Bulge, yet all W UMa-type systems discovered 
there appear to be Disk population objects distributed over the long line of sight to the Bulge 
QRucinski 1998b|) . While the W UMa-type binaries are very short-period systems, with orbital 



periods within 0.25 < P < 0.6 days (this class of binaries is normally assumed to have periods 
shorter than one day), use of a Fourier decomposition on all available OGLE-I light curves showed 
that main-sequence contact binaries continue to exist in the Galactic Disk population to orbital 
periods of about 1.3 — 1.5 days and then abruptly vanish, with very few contact systems having 



longer periods (Rucinski 1998a = Paper III, see Figure 2 in this paper). Only four apparently 
contact systems were discovered with the orbital periods P > 10 days. Because of the long periods, 
these could not be main sequence stars such as those in the typical W UMa-type systems, but 
giants. Since there exist no theoretical predictions nor any observational indications that contact 
systems can be formed with giant components, no firm explanation for the long-period systems in 
BW could be offered at that time of Paper III. However, it was noticed that their light variations 
could be rather easily be explained by orbital synchronous rotation of a single distorted component 
filling its Roche lobe. Thus, these long-period binaries could consist of just one evolved component 
and some other star which - because of its small size and/or low temperature - would be invisible, 
yet would exert a sufficiently strong gravitational pull to distort the visible component. 

In this paper we suggest that the "long-period contact binaries" are in fact representatives 
of the well known group of (5 Lyrae-type binaries. The primary basis for the argument is the 
well-populated sequence of long-period (LP) systems which has been established among contact 
binaries discovered by the OGLE-II project in the Small Magellanic Cloud (SMC)( [tJdalski et al. 



1998a). We have found members of this sequence by using the same Fourier light-curve shape 
filter as in the previous papers of this series. Although the filter is meant to select genuine contact 
binary systems, it can also pass binaries which are not in contact, but which show sufficiently 
strong ellipsoidal variations and shallow or absent eclipses. 

The hypothesis the LP systems are in fact binary systems with only one, Roche-lobe filling 
component dominating in the total light leads to a prediction of a period-luminosity-color (PLC) 
relation for such systems. Such a relation is expected to be very similar to what has been 



suggested for the W UMa-type stars ( Rucinski 1994 ), but with only one component defining the 



global relationship. If the total mass and the mass-ratio are not correlated with the orbital period 
(however, this is an assumption only), a PLC relation will simply result from the Kepler-law scaling 
of the total systemic dimensions (and thus of the stellar radius) with the orbital period, R oc P 2//3 . 
Further conversion to a luminosity dependence, L cx B?T^j, and then to an absolute- magnitude 



1 b.ttp: / /www. astrouw.edu.pl/~ogle/ or 



http: / /www. astro, princeton.edu/~ogle/ 
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dependence, will lead to My oc —10/3 log P. In this paper we show that a sample of (3 Lyrae 
systems in the SMC, selected with the exclusion of blue systems ((V — I) > 0.25), forms a 
surprisingly uniform and coherent group for which an absolute-magnitude relation - expressed as 
My = My (log P, color), with the expected logP dependence - can indeed be established. 

We define the sample in Section || and then show the period-color relation (Section ||) as the 
best tool to delineate the sequence of the long-period contact systems. In Section || we show that 
in fact these are (3 Lyrae- type binaries, not genuine contact systems. The color-magnitude diagram 
(Section |B|) leads directly to a period-luminosity-color which we express as a absolute-magnitude 
calibration (Section ||). We discuss the OGLE-II data for the systems discovered in the Baade's 
Window sample in Section |j] and give estimates of the spatial frequency of the LP systems in 
Section |8[ Section || summarizes the results of this paper. 



2. SELECTION OF THE SAMPLE 

The first step in the selection of the sample was finding the contact systems among all 
eclipsing binaries found by Udalski et al. 1998a| in the SMC. Then, following the results of the next 



Section [3|, which suggest a change in the sample characteristics at periods of about 10 days, the 
long-period representatives of the sequence were further selected. The 02 — 0,4 Fourier light-curve 



shape-filter has been applied to all eclipsing binaries listed in Udalski et al. 1998a. The operation 



of the filter is shown in Figure |l| while the details of the filter have been described before (Papers 



I, II, III and in Maceroni & Rucinski 199£ = Paper IV). The filter requires the cos(4</>) term to be, 
in its absolute value, relatively small in relation to the cos(2</>) term (4> is the orbital phase). A line 
on the ci2 — 04 relation, pre-computed for a case of the inner contact using light curve synthesis 
models, defines a locus of systems with components just touching. Contact systems show small 
(in absolute value) 04 coefficients when compared with 02- In effect, the Fourier light-curve filter 
selects eclipsing binaries with dominance of geometric distortion effects (which define the overall 
amplitudes and control the size of 02) over eclipse effects (contributing primarily to the 04 term). 

Normally, the Fourier filter works very well in correctly identifying contact binaries. It 
also rejects very effectively semi-detached binaries of the Algol type which show deep eclipses. 
However, the filter may give wrong results for a semi-detached binary with very shallow or absent 
eclipses, in the situation of a great disparity of sizes of the components, with a large visible star 
showing a considerable tidal distortionP]. Then the 04 coefficient can also be small relative to 02 
and the binary may be erroneously classified as a contact system. This is exactly what we suspect 
to be the case for the LP binaries described in this paper. 

As an aside, we note that the OGLE-II project eliminated small amplitude systems from 



2 The filter fails also for systems with secondary minima displaced from phase 0.5. Such systems had to be removed 
manually form the sample. 
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consideration at the early stage by exclusion of the so-called "ellipsoidal variables" (Udalski et al. 



1998a| ). This means that some small-amplitude contact, semi-detached and perhaps even detached 



binaries have been excluded at the stage of formation of the sample. We consider this choice as 
somewhat unfortunate because criteria for rejection were not clearly formulated in the discovery 
paper. Thus, the sample cannot be used for statistical studies involving light variation amplitudes 
or their derived quantities (such as ci2 and 04) at the (poorly defined) low-amplitude ends of their 
distributions. 



3. PERIOD COLOR RELATION 

The sequence of long-period contact systems is very well defined in the period-color (PC) 
relation which - as was many times stressed in this series - plays a special role for contact binary 
stars. The left panel of Figure |2| shows the PC relation for the all SMC systems selected as contact 
with the 02 — 04 light-curve shape filter. The sequence consists of 38 systems. It starts at 11 days 
(in fact a couple of systems with periods of 8 - 10 days may also belong to the sequence) and 
continues to 181 days. It has an inverted slope relative to the one normally observed for most 
common contact binaries of the W UMa-type: Typical contact binaries become hotter as their 
periods get longer, a trend which is explained by a progression upward along the main sequence. 
The inverted sense of the slope for the LP sequence suggests that they become cooler with the 
increase of dimensions, a property normally associated with giants. 

The SMC sample is not the first one to show the LP sequence: In addition to the indications 
from the OGLE-I Baade's Window (BW) sample that we mentioned in Section |], the MACHO 



project has also detected several similar systems in LMC; the latter have been analyzed in Rucinski 



1999| . We show the LMC systems in the right panel of Figure y. Although the BW and the LMC 
samples give a qualitative support to what will be discussed in this paper, we cannot utilize the 
BW and LMC data together with the SMC data for the following reasons: (1) It is has been found 
that the OGLE-I photometry for BW suffers from zero-point offsets as well as from scale errors 
resulting from the non-linearity of the CCD detector. We discuss the details and present limited 
new OGLE-II data for three systems in Section [?| (2) The MACHO sample used unconventional 
spectral bandpasses which resulted in photometry only approximately transformable into the 
(V, V — R) color system. Thus, neither the OGLE-I BW sample nor the MACHO sample can be 
compared with the new, high quality OGLE-II data for SMC. 



4. LIGHT CURVES 



The OGLE-II data for the SMC variable stars consist of light curves in the / band and of the 
maximum light magnitudes and color indices {B — V) and (V — I). The light curves are shown 
in Figures || and ^| while the photometric data are given in Table |l[ The first columns of the 
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table contain the OGLE project data: First, the identification numbers (the consecutive ones used 
in this paper and then the OGLE-II numbers), followed by the orbital periods in days and the 
three photometric quantities, the I ma x magnitude and the color indices (B — V) and (V — I). The 
OGLE-II identification number is made of the SMC field number (integer part of the number) and 
of the star number (six digits after the period). Among the 38 systems found in the SMC fields, 
three have additional identifications in the neighboring fields. The collection of the light curves in 
Figures [3| - 1| shows only the light-curve data based on the primary identifications. The light-curve 
Fourier decomposition coefficients a,, which are given in the last six columns of Table |l|, have 
been calculated for both identifications, but we used the primary ones in selection of the contact 
binary systems through the 02 — 04 Fourier filter. As in the previous papers of this series, a% are 
the cosine coefficients while b\ is the only sine term representing light curve asymmetries. 

The light curves shown in Figures || and [| are not of typical contact systems with very similar 
depths of eclipses. All show unequally deep eclipses which are characteristic for (3 Lyrae (EB) light 
curves. The light curves are well defined for most systems. Twenty seven systems with relatively 
large amplitudes of light variations have been selected to enter with full weight in determinations 
of luminosty calibrations in Section ^. Eleven systems with poorer light curves are marked by 
asterisks in Figures || - || and in Table ^ they are marked by open circles in the remaining figures. 

The similarity of all LP systems in terms of their light curves is confirmed by the a<i — a\ plot in 
Figure |5] where we see a relatively tight relation, almost a proportionality between the coefficients. 
For normal W UMa-type systems, the a\ coefficients are confined to a band —0.02 < a\ < 0, 
without any clear correlation with ai (see Rucinski 1997b| , Figure 5 there). Since the coefficients 



are defined in terms of light intensity units, the proportionality of a\ to a<i implies similar light 
curves in magnitudes, with similar differences of the eclipse depths. This is exactly the property 
that was noticed when an attempt was made to model the light curve of the LP system in the BW 
sample, BW0.036 (Paper III): For a system where the only visible component fills the Roche lobe, 
the relative depths of eclipses are expected to be very similar almost irrespective of the orbital 
inclination and only very slightly depend on the mass-ratio. 

Although the light curves are similar, they also show some small differences in that some 
show narrow secondary eclipses while some show relatively more scatter at both maxima. In 
one case, LP#9, we see well defined eclipses indicating that both components contribute to the 
systemic light; we show later that this is exactly the system which deviates from the PLC relation. 
We have made no attempt to model the light curves of the LP systems. We feel that without 
spectroscopic mass ratios, this would be futile exercise. We note that color curves cannot be 
currently formed because of much smaller numbers of the B and V observations relative to the 
number of I observations. However, all the light curve features can be explained by the same 
model which is applicable to the well-known prototype of (3 Lyrae. 



At this point, we would like to make a clear distinction between the light-curve shape 
designated by EB from the name of close-binary stars [3 Lyrae. From now on, (3 Lyrae will be used 
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here as a binary-system prototype, not as a designator of of the EB light-curve shape. That is, 
when calling a star a (3 Lyrae type, we will consider a binary with only one, giant-type component 
which dominates the total light, the other component being a main sequence star or collapsed 
star (perhaps surrounded by an accretion disk) with only a minor influence on the total systemic 
brightness and its light curve. Such a model is now well established for the prototype, following 



the extensive research of /3 Lyrae and of the related W Serpentis-type binaries ( Huang 1963 , Woolf 
1965j IWilson 197| , |Plavec 19851 , |Hubeny fc Plavec 1991| , |Harmanec fc Scholtz 1991 , pe Greve k 



Linnell 1994 ) and of various X-ray sources with giant-star mass donors. We stress that in calling a 



star a (3 Lyrae binary, we focus on two essential properties: the semi-detached model and the giant 
characteristics of the Roche- lobe filling component. Later on, to establish an absolute- magnitude 
calibration, we add the third property: the sufficiently yellow/red color which, in most cases, 
indicates a negligible photometric importance of the other component. 

The light curves of the LP (3 Lyrae systems in the SMC sample appear to be symmetric, 
with identical heights of the light maxima (i.e. the so-called O'Connell effect is nonexistent for 
these systems). We also see no correlation between the sine coefficients, b\, and the coefficients 
measuring the eclipse-depth differences, a\. Such a correlation, manifesting itself as a preference for 
higher first maxima with the increased eclipse-depth difference, was observed for the short-period 



EB systems (P < 1 day), appearing in the Baade's Window sample and analyzed in Rucinski 



1997b| . An explanation offered there was that these are semi-detached binaries with matter being 



transferred from the more massive, hotter to the less-massive, cooler components. Such systems 
are very different from the systems discussed in this paper in that they do not have to contain 
a giant, Roche-lobe filling component. In fact, judging by short orbital periods, many of them 
consist of main sequence stars first time evolving into the mass exchange phase or, perhaps, 
experiencing the semi-detached stage of the thermal-cycle relaxation oscillations ( Lucy 1976| , 



Flannery 1976, Robertson k Eggleton 1977, Lucy k Wilson 1979). 



Finally, we note that the light curve amplitudes of the LP systems, as measured by the a,2 
Fourier coefficient, do not show any correlation with the orbital period or the color index of the 
system. Note also that the light curve shapes in Figures || - || do not correlate in any obvious way 
with the periods (the periods are given in parentheses above each panel in these figures). In this 
respect, the sample appears to be uniform in the sense of consisting of very similar systems. 



COLOR MAGNITUDE DIAGRAM 



The SMC data are an excellent material for forming a color-magnitude diagram (CMD) 
and for studying location of the LP systems in relation to the well recognized stages of stellar 
evolution. The distance modulus is relatively well known and will be soon known even better 
once the matter of systematic deviations between various methods of distance determinations to 



Magellanic Clouds is resolved QUdalski 19981 |Udalski et al. 1998b| , pibson 1999| , [Nelson et al 



2000j |Popowski 2000|) . Also, the average reddening for the SMC is small. In contrast, the Baade's 
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Window sample (Paper III) suffered from large and uncertain reddening corrections and from 
the (then unrecognized) photometric errors, more fully described in Section [7|. In addition, the 
distances of the BW systems were of course different precluding formation of an independent 
color-magnitude diagram (other than the one involving an absolute-magnitude calibration). In 
principle, the MACHO sample for the LMC ( [Rucinski 1999 ) could be used to form a CMD, but 



the data were in a non-standard photometric system and the reddening corrections (and their 
variations) are larger than typically for the SMC. 

Figure ^ gives the observational CMD for the LP Lyrae systems in the SMC, with the 
main stellar density enhancements (such as the main sequence or the red clump) marked by 
shading. The approximate discovery completeness limits of 90% (the slanting line starting at 



/ = 17.5 and V — I = 0) and of close to 0% (1.5 mag. below), following Udalski et al. 1998a, are 
shown by slanting, dotted lines. In transforming from the observed, maximum-light photometric 
quantities to the de-reddened, absolute magnitudes, we used: (m — M)q = 18.8, Ey-i = 0.12 and 
Ay = 2.5 Ev-i- The LP systems are much more luminous than the main-sequence stars. They 
are also systematically brighter than typical giants in the SMC; this may be partly a selection 
effect as brightest stars are expected to have better photometry. Comparison with in evolutionary 



tracks in pVlassey et al. 1995| indicates that the LP binaries considered in this paper cover the 
area normally associated with evolved, moderately massive (3 — 5Mq) stars within the range of 
effective temperatures of 3.60 < logT e jf < 4.00. Later on, we set a blue limit at (V — I)o = 0.25 
so that the effective temperature range is cut from above at logT e ff ~ 3.9. 

The orbital periods of the LP j3 Lyrae systems are shown in Figure | by the size of the 
symbols. One immediately sees a period progression in that the systems with the longest periods 
appear in the upper part of the CMD. The longest periods, P > 100 days, appear exclusively 
among red giants above the Red Clump location, but a similar period progression is also visible 
in the upper-left, blue part of Figure |6] where five systems of similar color index have luminosities 
scaling with their orbital periods. This type of the period dependence indicates that a period - 
luminosity - color relation can be established for the LP systems. 



6. PERIOD LUMINOSITY COLOR RELATION 

6.1. Calibration of the relation 

The previously established absolute magnitude calibrations for contact binary systems of 
the W UMa-type are of limited value here, primarily because of the vastly different ranges of 
the orbital periods. While the W UMa calibrations are applicable to the range 0.25 - 0.6 days 
( Rucinski Duerbeck 1997| )P|, the calibration that we would like to establish will cover a large 



3 A blind application of the (V — I)-based calibration of 
My estimates about 3.5 magnitudes too bright. 



Rucinski & Duerbeck 1997 



to the SMC LP systems gives 
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range of the orbital periods from 11 days to 181 days. 

We set aside and do not discuss the contact systems with the intermediate periods, starting 
around 0.6 - 1.5 days at the short end and extending to 10 days. Such systems apparently 



do not exists in the solar neighborhood and are exceedingly rare in the BW sample (Rucinski 



1998a ), but - possibly due to the discovery selection effect preferring brightest objects - are quite 
common in the present SMC and in the MACHO LMC samples (Figure |J). We feel that the 
intermediate-period (1.5 — 10 days) systems which are visible in the SMC and LMC samples are a 
mix of genuine contact systems with the moderately-short representatives of the (5 Lyrae systems 
discussed in this paper. They will be the subject of a separate investigation. 

Assuming the same values for the distance modulus and reddening as in the previous section, 
the linear regression for the whole available data set of 38 systems (with 11 systems given half 
weights) gives: 

My = -3.90 log P + 3.27 (V - J) + 2.40, a = 0.50 (1) 
My = -4.05 log P + 3.01 (B - V) + 3.46, a = 0.63 (2) 

The standard errors (per data point) a are large and one may question the validity of linear fits 
in so poorly constrained cases. The deviations from the linear regressions are shown in the upper 
two panels of Figure 0. Note, that the bluest stars appear to deviate the most from the linear 
dependencies. If we exclude these systems setting the limit at (V — 1)0 > 0.25, then the linear fits 
improve substantially reaching the level of precision of about 0.35 mag. This level of accuracy is 
similar to that for W UMa-type binaries and for Cepheid variables. For the limited data set of 33 
yellow-red systems, the linear weighted regression equations are: 

My = -3.43 log P + 2.04 (V-I) + 2.80, a = 0.34 (3) 
My = -3.30 log P + 1.49 (B-V) + 3.38, a = 0.38 (4) 

The standard deviations are now quite respectably small so that the equations can be used for 
prediction of absolute magnitudes. The deviations from the linear fits for this case are shown in 
Figure ||. Realistic uncertainties of the coefficients for the sample of 33 yellow-red systems have 
been established by application of a bootstrap re-sampling process. The median values and the 
1— a ranges give the following relations: 

My = -3.411°;^ logP + 2.00tg|?(y- J) + 2.80l°;P (5) 
My = -3.29lg;itlogP + 1.48l8:ii(B-VOo + 3.36l8:g (6) 



The log P coefficients are very close what we suggested would be the case on the basis of 
very simplified considerations in Section [l], My oc —3.33 logP, which is highly gratifying. But we 
achieved that by removing the five blue systems with (V — 1)0 < 0.25. This removal is not entirely 
arbitrary; we discuss it in the next subsection. One can fix the period coefficient at the expected 
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Kepler-law value of —3.33 and determine the other two regression coefficients with higher much 
accuracy. The results for 33 yellow-red systems are: 



My + 3.33 log P = 1.94 (V — I)o + 2.72, a = 0.34 (7) 
My + 3.33 log P = 1.51 (B — V)o + 3.40, a = 0.38 (8) 

while the bootstrap-estimated median values and the 1—a ranges lead to: 

My + 3.33 log P = 1.94^ (V - J) + 2.72±g;J| (9) 
My + 3.33 log P = 1.501°;^ (B - V) + 3.41^ (10) 

We suggest that until better and more extensive data become available, Eqs. (^)-@ or 
Eqs. (|7|)-(||) serve as the preliminary calibrations of the PLC-/? relation. 



6.2. Limitations of the PLC-/3 relation 

One of the limitations of the PLC-/3 relation is the necessity to limit selection of the systems 
to redder than (V — J)o = 0.25. There a few reasons why the blue stars may deviate from the 
general trends: 

• The underpinning of the PLC-/? relation is the assumption that the large, Roche-lobe filling 
component dominates the photometric characteristics of the system. Only then one can 
expect a tight relation between the size of this component (governed by the Roche lobe 
scaled with the size of the whole system) and the total luminosity of the binary system. 
Even a relatively small, but hot and hence possibly luminous companion can modify the 
color index to make it blue. Thus, blue color index may be considered as a direct indication 
that the hot companion of the Roche-filling star is contributing substantially to the system 
light. The well-observed system LP#9 (OGLE-11^3. 134555) which shows clear eclipses is an 
excellent illustration of the point in that it is one of the blue systems removed from the final 
calibration. 

• It is well known that the optical-range color indices, such as (B — V) or (V — I), lose 
sensitivity to effective temperature for very hot stars. The serious problems in the color 



- effective temperature transformations occur at logT e ff > 4.5 (Massey et al. 1995). In 
our case, the blue stars which deviate from the PLC-/3 relation are not that hot, and have 
temperatures within the range 3.9 < logT e jj < 4.0. 

P Lyrae is itself a good case of a blue binary system which does not obey our PLC-/? 
relation. In this case, it is not the contribution of the other component to the total light, 



but apparently the Roche-lobe filling component is too luminous. Dobias & Plavec 1985 
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summarize the parameters of relevance to our discussion: P = 12.937 days, My = —4.7 for 
both components and My s = —4.1 for the visible component, {B — V)q = —0.13 (which 
implies (V — I)o ~ —0.1. Our calibrations predict My 1 = —1.3 and My V = —0.5 for the 
visible component. This component is thus some 2.8 - 3.6 mag. over-luminous relative to 
the calibrations. 

• Possibly, use of single linear relations for the color dependence is simply inappropriate over 
the whole wide range < (V — J)q < 1.6 observed for the LP (5 Lyrae systems. The five hot 
systems in our sample do show the steep relation between My and log P. If we lump all five 
and use their average color indices ((V — I)o = +0.109 and (B — V)q = —0.067) to enter 
these in Eqs. @ and (||), then My oc —3.10 logP, with the zero points, +0.87 and +1.19 for 
the (V — I) and {B — V) based relations, respectively. These zero points are very different 
than for the more common yellow-red LP systems, leading to relations by about 2 to 2.5 
magnitudes brighter. This would go some way towards expaining the over-luminosity of the 
(3 Lyrae less-massive component. 

A serious limitation of any approach requiring de-reddened colour indices is the uncertainty 
of the reddening correction. Although the colour-index terms have relatively small coefficients 
and the strongest dependence is in the logP term, nevertheless the reddening corrections must 
be applied to the observed data. For the SMC OGLE-II data, two indices are available, (B — V) 
and (V — I). Can one use form a reddening- free index, similar to the well-known Q index? 
Unfortunately, this is not the case for the particular combination of the indices. An index 
constructed in the form: Q' = (B — V) — k (V — I), where k = Eb-v / 'Ev-i is n °t a monotonic 
function of the indices: It shows a shallow increase from (V — I)q = —0.2 to 0.9, then a broad 
maximum around 0.9 < (V — /)o < 1-4, and a steep decline further to the red. Thus, while the 
conventional Q-index technique based on the (U — B) and (B — V) would work, the one based on 
the (B — V) and (V — I) indices would fail. The reddening correction uncertainty will remain the 
main limitation of the PLC-/3 relation. 



6.3. External checks: The LMC and M31 data 

Undoubtedly, calibrations of the PLC-/? relation will be substantially improved when similar 
OGLE-II data for the Large Magellanic Cloud become available. At present, support from other 
sources is weak. As has been pointed before, the currently available MACHO data for the LP 
systems in LMC (|Rucinski 1999|) are in a non-standard photometric systems only approximating 



the VR filter system. Seven systems in this sample have orbital periods longer than 10 days, 
in the range 12 < P < 52 days. Unfortunately, all seven are quite blue (V — R < 0.35, 
see their location in right panel of Figure ^). It is nevertheless possible to establish a tight 
linear relation, similar to the calibration Equations @ - (10), for these systems. The relation 



(My = -0.26 logP + 4.31 (V - R) - 1.58 for m - M = 18.3 and Ey_ R = 0.2) has a small 
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dispersion a = 0.1 mag., but we do not believe that this relation is of much validity because the 
photometric system is non-standard and all systems are blue. 

The data for binary systems in the Andromeda galaxy detected by the DIRECT project also 
do not provide an external check on the calibration of the PLC-/3 relation. The project found four 
LP systems of f3 Lyrae-type in M31 ( [Stanek et al. 1998 , Stanek et al. 1999. Kaluzny et al. 1999). 



All four have periods within 11 to 14 days, so that they do not probe the whole span of the orbital 
periods of the SMC sample. The DIRECT systems are: D31A-6423, D31C-13944, D31D-5186, 
D31D-7221 (we will call them below by consecutive numbers, 1-4). The systems $4 and #4 
have very poor photometric data (especially the crucial color indices are poorly determined) while 
the system #2 has too blue color index (uncorrected V — I = 0.13) to belong to our proper LP 
sequence. Only the third system, D31D-5186 would qualify, with P = 11.80 days, V max = 19.5, 
V — I = 0.41 and the very uncertain B — V = —0.2. The directly evaluated value of My, assuming 
m — M = 24.5 ( Stanek Garnavich 199§| ) and Eb~v = 04 is My = —5.3, while the one estimated 



from Eq. (||), assuming the same reddening, is My = —0.2. Thus, the discrepancy between the 
two estimates is very large, either shedding doubt on our SMC calibration or indicating that the 
companion is over-luminous (as in the very case of (5 Lyrae). We note that the DIRECT project 
can discover only the very brightest binaries in M31 and these would be the ones which would be 
most likely the ones which deviate the most from the calibration. Also, the color used is highly 
dependent on the assumed reddening which is only guessed in this case. We cannot exclude a 
possibility that (V — I)o for the system ^3 is smaller than our threshold of 0.25; then the system 
would belong to the blue group which we explicitly excluded from consideration when establishing 
the PLC-/3 relation. 



7. LONG-PERIOD SYSTEMS IN BAADE'S WINDOW 

It has been realized recently that the OGLE-I photometry for Baade's Window (BW) suffered 
from systematic zero point offsets and linearity problems resulting from the non-linearity of 
the original CCD detector. At the brightness level of the Bulge Red Clump Giants, the offsets 



are moderate: 0.021 in V, 0.035 in / and thus 0.056 in (V — I) (Paczynski et al. 1999). More 
worrisome are however the non-linearities, relatively larger in the ^-filter data for red stars 
and affecting particularly strongly the (V — I) color indices of faint and red stars. Dr. Udalski 
(private communication) estimated that deviations in (V — I) may be as large as 0.2 mag. for 
faintest surveyed stars. Since the color indices play an important role in our absolute-magnitude 
calibrations, the deficiency of the OGLE-I photometry will require a more extensive re-discussion 
of the results of Papers II and III. 

The OGLE project, currently in its Phase II, is characterizing the systematic errors of 
OGLE-I. Some stars already have excellent new data in V, {B — V) and (V — I). At our request, 
Dr. Udalski kindly checked the availability of new data for the sequence of the systems with 
P > 1 days which had been discussed in Rucinski 1998a| (Paper III). About two thirds of these 
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stars have new photometry. Among the four BW systems with P > 10 days, three have new 
photometric data. These have been made available to us in terms of the average light levels, 
requiring additional corrections to the maximum light level. They are based on typically 30 
observations in / and 10 observations in V. We applied the differences AI = I\ — 1%, as found on 
the basis of the OGLE-I light curves (Paper III), to convert the data from the average (I2) to the 
maximum (I2) light levels. The details are given in Table ||. 

The left panel of Figure |2| shows the OGLE-II data for the three LP (3 Lyrae systems in the 
period-color plane, with the assumed reddening corrections the same as in Paper III (the three 
systems are marked by filled asterisks). The corrections were found on the basis of an assumption 
that the reddened stars are further than 2 kpc and thus suffer the same amount of reddening 
as the Bulge stars in the BW ( |Stanek 1996 ). This assumption can be checked for consistency 
using the calibration in Eq. (^). The estimates of the distances are given in the last column of 
Table ||. The errors reflect the uncertainty of the calibration equation of about 0.3 mag. The 
results are consistent with two systems in the Bulge (possibly on the other side of it) and one 
system possibly in front of the Bulge or inside its elongated part directed toward us. The fourth 
system, #0.036, without the new data, also appears to be located in the Bulge. We treat these 
distance determinations as very preliminary and disregard a possible metallicity mismatch between 
the calibration established for the SMC systems and the possibly metal-rich systems in Baade's 
Window. 



8. SPATIAL FREQUENCY 

The spatial frequency of the long-period (3 Lyrae-type systems is difficult to estimate because 
it is practically impossible to estimate the OGLE-II search volume in the SMC. In terms of 
the intrinsic luminosities, the LP (3 Lyrae-type systems occupy the range —5 < My < +1 (the 
calibrated part extends to My ~ —3) so that they can be compared with the Cepheid variables. 
The OGLE-II project observed 2155 Cepheids in the SMC (Dr. A. Udalski - private information) 
which can be compared with 38 long-period (3 Lyrae-type systems discussed here. While a large 
fraction were previous discoveries and the LP binary systems are new, this gives a rough idea 
about the frequency of occurrence of the LP (3 Lyrae-type more than 50 times below the spatial 
frequency of the Cepheids. As an aside, notice when comparing the two types of variables, that 
the period term in the LPC-/3 relation appears to be close to the expected My oc —3.33 log P, 



while the Cepheids show a period dependence, with the period coefficient —2.78 ( Feast fc Walker 



1987 ) or —2.76 ( Madore Sz Friedman 1991 ). The similarity of these dependencies should not be 
surprising as both are driven by the same dimension - response time relationship of the dynamical 
time scale. 

The Baade's Window data give us only a very approximate constraint on their frequency 
in our Galaxy. If one determines the average frequency of occurrence by simply dividing their 
number (four) by the volume of the OGLE-I search in the BW area (about 2.5 x 10 7 pc 3 to the 
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Galactic center), then an estimate of the frequency is 1.6 x 10~ 7 LP systems per cubic parsec. 
However, the LP /? Lyrae systems appear to be strongly concentrated towards the Galactic Bulge 
so that such a naive estimate does not take into account the increase in numbers of stars as the 
line of sight penetrates the Galactic Bulge. We recall that the W UMa-type systems of the disk 
population approximately follow the stellar density in the BW direction with the average density 
of about 7.6 x 10~ 5 systems per cubic parsec dRucinski 1998b|) . 



9. DISCUSSION 

This paper discusses consequence and observational support for an assumption that the 
long-period (P > 10 days), apparently contact binary systems observed in Baade's Window, 
and in LMC and SMC, are in fact (3 Lyrae systems with the cool components dominating their 
total luminosity. We analyzed in detail the LP systems in the SMC and found that they form 
a well-defined, moderately tight, "inverted" period-color relation, implying a color-size relation 
characteristic for giants rather than for main-sequence stars. The principal support, however, 
comes from the existence of a period-color-luminosity relation for those which are sufficiently cool. 
Tentatively, we set the color-index limit at (V — I)o > 0.25. An absolute-magnitude calibration 
of the new PLC-/3 relation in the form My = My (log P, color), established in a preliminary form 
on the basis of 33 systems in the SMC, has the expected slope My cx —3.33 log P. It offers an 
exciting application for distance determination to nearby galaxies. The LP (3 Lyrae systems are 
the brightest among binary stars and their continuous light variations make them also very easy to 
spot. They occupy roughly the same range of the absolute magnitudes as the Cepheid variables, 
but appear to be by about 50 times less common in space. Therefore, the PLC-/3 relation can play 
only an auxiliary - although entirely independent - role in distance determinations. 

We see sequences of ostensibly contact, long-period systems not only in the SMC, but also in 
Baade's Window and in MACHO data for LMC. However, the latter data sets could not be used to 
improve the absolute magnitude calibration of the PLC-/3 relation for reasons described in detail 
in Sections |3| and [7[ It is expected that the data for LMC currently collected by the OGLE-II 
project will provide an excellent check and improvement of the tentative results presented here. 

We conclude by observing that the j3 Lyrae binaries discussed in this paper are expected 
to show some similarities to the cataclysmic variable (CV) binaries. The CV binaries are also 
known to have secondary (less massive components) filling their Roche lobes. As was observed by 



Faulkner et al. 1972 , a combination of Kepler's law, a 3 oc (M\ + M2) P 2 , with the approximate 
relation for the relative radius of the Roche lobe filling star, P2/0 = 0.462(^/(1 + g)) 1 / 3 , with 
q = M2/M1 < 1, gives a relation similar to that for pulsating stars: P = 0A38/^/~p, where the 
period P is in days and the mean density of the star, p, is in g cm~ 3 . This relation is valid for 
small mass ratios (q < 0.5), but corrections for larger values of q are small. Because the range of 
the observed period for the LP sequence in the SMC is large, one expects the Roche-lobe filling 
components to have mean densities within a large range from about 2 x 10 -3 g cm -3 to less than 
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10 5 g cm 3 , that is within typical range observed for yellow-red giants and supergiants. 

Special thanks are due to Dr. Andrzej Udalski for his prompt response - in spite of many 
ongoing activities - to a request to extract new OGLE-II data for the LP systems observed by 
OGLE-I project in Baade's Window and for several explanatory comments on the data. Thanks 
are also due to Dr. Stefan Mochnacki for very useful comments and suggestions. The authors 
would like to thank the OGLE team for making the data freely available over the computer 
networks. 
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Captions to figures: 

Fig. 1. — The Fourier light-curve shape filter selects contact binaries by comparison of the a-i and 
a4 cosine coefficients. Systems with properties of contact binaries are located below the curved line. 
The small dots give the data for all eclipsing systems in the SMC while the large filled circles mark 
38 long-period (P > 10 days) contact systems discussed in this paper. The open circles are used 
for 11 systems with relatively poorer light curves which have been given half weight in luminosity 
calibrations discussed further in the paper. Note that contrary to the results for Baade's Window 
(Rucinski 1997a), the contact domain is rather thinly populated in the SMC, mostly due to the 



absence of the genuine W UMa-type binaries (P < 1 day). We give arguments in the paper that 
the apparently contact long-period systems are in fact semi-detached binaries of the j3 Lyrae-type 
with component-distortion effects dominating over shallow or absent eclipses. 

Fig. 2. — The period-color diagram for all contact systems in the SMC (this paper, left panel) 
and LMC (the MACHO data, right panel, as discussed in Rucinski 1999). The large circles show 
systems discussed in this paper; the open circles mark those which have been assigned half weights 
in luminosity calibrations discussed later on. The SMC data are based on the (V — I) color index 
while the LMC data are based on the (V — R) color index. The curves give the blue short-period 
envelope for normal, short-period W UMa-type systems. For simplicity, we have assumed the same 
reddening vectors for all systems in the SMC, Ey-i = 0.12 (vertical lines). In subsequent sections 
we find that systems with (V — 1)0 > 0.25 provide a well defined period - luminosity - color 
relation. Three large asterisks show locations of the long-period systems in Baade's Window; they 
are more fully discussed in Section [7[ We cannot directly utilize the BW or LMC data in this 
paper, but point out the existence of the LP systems which appear to be in contact and form a 
sequence extending to long orbital periods. 

Fig. 3. — A collection of / band light curves for 38 long period contact systems in the SMC. 
This figure contains 24 light curves; the remaining curves are in the next figure. The systems are 
identified by numbers above each of the panels with the orbital periods given in parentheses. An 
asterisk added to the star number indicates half-weight given in the luminosity calibrations. The 
vertical scale of each panel spans one magnitude and extends from 0.6 mag. below to 0.4 mag. 
above the average magnitude. The maximum light magnitudes, the orbital periods, the color data 
as well as the Fourier cosine coefficients are given in Table 1. 

Fig. 4. — The same as for Figure 3, but for 14 among 38 long period contact systems in the SMC. 

Fig. 5. — The first cosine term of the light-curve Fourier decomposition, a%, appears to scale in 
proportion to the largest term, 02- This confirms that the light curves are very similar in shape. 
The open circles mark systems with half weights in the calibration solutions. 

Fig. 6. — The color-magnitude diagram for the LP systems in SMC. The orbital periods of 
individual systems are marked by the size of the circle. The right and upper axes show the de- 
reddened color index and absolute magnitude data for Ey-i = 0.12 and (m — M)o = 18.8. The 
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slanted dotted lines give the approximate locations of the 90% and 0% completeness levels of the 
OGLE-II variable star search. The broken lines represent the constant period lines according to 
the calibration of our PLC-/3 relation as in Eq. (7). They extend to the red of (V — I)o = 0.25. The 
current calibrations are not valid to the blue of that point (marked by the vertical dotted line), 
although a period dependence with a correct slope is observed for such systems as well. The main 
sequence (MS) and the red clump (RC) of the SMC are marked by shading. 

Fig. 7. — The deviations from the regression lines of the calibration in Equations (1) - (2). Note 
the large negative deviations at blue color indices. The labels are the consecutive numbers in 
Table 1; they are given for systems deviating more than 0.75 mag. from the linear regression. Note 
in particular the location of the system LP#9 which is discussed in the text. The dotted lines give 
the color-index cutoff lines at (V — I)q = 0.25 and the corresponding (B — V)q = 0.21. 

Fig. 8. — Same as in Figure 7, but with a limitation to systems with (V — 1)0 > 0.25, as in 
Equations (3) - (4). Notice a small inconsistency - which we disregarded in our determinations of 
the PLC calibrations - in that 3 systems are slightly bluer than (B — V)q = 0.21 which roughly 
corresponds to the cutoff at (V — I)q = 0.25. 
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Table 1. [3 Lyrae systems with P > 10 days in the Small Magellanic Cloud 





LP# 


SMC# 


P(d) 


^max 


(B-V) 


(V-I) 


ao 


"i 


a 2 


0,3 


0-4, 


bi 




1 


1.012977 


60.1884 


14.279 


0.020 


0.257 


0.859 


-0.043 


-0.129 


-0.024 


-0.011 


-0.005 




2 


1.035641 


29.1173 


17.711 


0.378 


0.735 


0.936 


-0.015 


-0.068 


-0.002 


0.002 


-0.002 




3* 


1.062436 


13.6637 


18.351 


0.768 


1.039 


0.972 


-0.022 


-0.041 


0.003 


0.004 


-0.007 




4 


3.002900 


89.8238 


16.183 


1.251 


1.372 


0.939 


-0.011 


-0.064 


-0.011 


0.002 


-0.006 




5 


3.060633 


28.2032 


17.894 


0.752 


0.988 


0.873 


-0.025 


-0.114 


-0.015 


-0.022 


-0.004 




6* 


3.075910 


72.9871 


16.928 


1.003 


1.177 


0.956 


-0.032 


-0.029 


0.000 


0.000 


-0.003 


o : 


7 


3.116114 


24.3557 


17.953 


0.386 


0.795 


0.921 


-0.019 


-0.089 


-0.017 


-0.008 


-0.002 


o . 


8* 


3.116383 


13.2463 


19.070 


0.216 


0.827 


0.968 


0.018 


-0.082 


0.006 


-0.011 


-0.003 


o ■ 


9 


3.134555 


11.8497 


16.283 


0.129 


0.333 


0.845 


-0.072 


-0.173 


-0.028 


-0.039 


-0.005 


(N ' 


10* 


3.134806 


55.5766 


17.898 




1.310 


0.899 


-0.019 


-0.099 


-0.011 


-0.011 


-0.001 


G ■ 


11* 


3.149100 


33.9911 


18.028 


0.750 


1.023 


0.937 


-0.003 


-0.050 


0.000 


-0.002 


0.008 


3 ; 


12 


3.193411 


28.4187 


17.197 


0.195 


0.555 


0.917 


-0.028 


-0.112 


-0.016 


-0.026 


0.003 


hj ; 


13 


3.238378 


43.8842 


14.880 


-0.047 


0.158 


0.957 


-0.009 


-0.053 


-0.008 


-0.001 


-0.005 


\o ■ 


14 


4.156495 


22.2481 


17.174 


0.219 


0.435 


0.899 


-0.025 


-0.123 


-0.011 


-0.023 


0.001 




15* 


4.179449 


20.3446 


18.088 


0.987 


1.275 


0.951 


-0.020 


-0.055 


-0.006 


-0.005 


-0.005 




16 


5.145588 


138.5914 


15.730 


1.056 


1.349 


0.881 


-0.044 


-0.122 


-0.028 


-0.025 


0.000 


i— i ■ 


17 


5.162500 


181.0044 


15.729 


1.583 


1.638 


0.928 


-0.021 


-0.083 


-0.022 


-0.010 


-0.002 


> ■ 


18 


5.180039 


97.1930 


16.065 


0.769 


1.041 


0.907 


-0.017 


-0.074 


-0.006 


-0.007 


0.000 


vo : 


19* 


5.257236 


11.7916 


18.788 


0.347 


0.778 


0.915 


-0.016 


-0.110 


-0.018 


-0.019 


-0.009 


l/000623 




(6.006812) 


11.7957 


18.790 


0.346 


0.755 


0.934 


-0.027 


-0.074 


-0.014 


-0.005 


-0.010 


20 


6.018257 


15.5577 


18.609 


0.536 


0.805 


0.883 


-0.035 


-0.150 


-0.023 


-0.032 


-0.003 


21 


6.094503 


112.4580 


15.290 


0.799 


0.922 


0.893 


-0.040 


-0.113 


-0.031 


-0.021 


0.009 


22 


6.110581 


88.3550 


16.046 


1.067 


1.171 


0.871 


-0.056 


-0.162 


-0.036 


-0.027 


-0.005 


23 


6.122285 


28.5879 


15.899 


0.008 


0.117 


0.905 


-0.012 


-0.095 


-0.009 


-0.005 


0.000 


24* 


6.199624 


62.1454 


16.547 


1.000 


1.067 


0.972 


-0.019 


-0.030 


0.001 


-0.001 


0.002 


Or 


25 


6.232915 


13.1408 


18.136 


0.338 


0.422 


0.897 


-0.032 


-0.118 


-0.014 


-0.017 


0.006 


26 


6.292000 


21.8606 


18.049 


0.916 


1.147 


0.937 


-0.049 


-0.104 


-0.003 


-0.020 


-0.019 


6 ■ 


27 


6.297249 


21.1419 


18.082 


0.406 


0.850 


0.923 


-0.023 


-0.089 


-0.014 


-0.004 


-0.005 


a: 


28 


7.018931 


12.6850 


18.497 


0.352 


0.599 


0.905 


-0.023 


-0.110 


-0.023 


-0.019 


-0.001 




29 


7.255665 


47.3855 


16.626 


0.813 


1.066 


0.875 


-0.037 


-0.156 


-0.027 


-0.035 


-0.003 


ctf ■ 




(8.030691) 


47.4408 


16.628 




1.053 


0.861 


-0.026 


-0.152 


-0.026 


-0.038 


-0.008 


> : 


30 


7.263170 


35.6454 


17.343 


0.576 


0.902 


0.858 


-0.042 


-0.156 


-0.021 


-0.040 


0.004 






(8.038913) 


35.6472 


17.352 


0.583 


0.884 


0.865 


-0.032 


-0.162 


-0.011 


-0.048 


0.006 




31* 


8.000105 


31.8424 


17.068 


0.817 


1.075 


0.970 


-0.012 


-0.042 


-0.004 


-0.004 


-0.012 


:ar 


32 


8.015451 


29.4313 


18.226 


0.737 


0.796 


0.873 


-0.033 


-0.152 


-0.016 


-0.042 


0.003 


33 


9.059110 


16.3410 


17.739 


0.657 


0.905 


0.924 


-0.012 


-0.094 


-0.010 


-0.012 


-0.005 




34* 


9.061975 


25.5624 


18.100 


0.760 


1.066 


0.981 


-0.037 


-0.042 


0.002 


-0.007 


-0.004 




35 


10.059436 


52.8010 


17.313 


0.946 


1.182 


0.906 


-0.010 


-0.116 


-0.004 


-0.026 


-0.002 




36 


10.078655 


129.0291 


15.308 


0.980 


1.156 


0.881 


-0.050 


-0.141 


-0.029 


-0.031 


-0.002 




37 


10.088122 


14.3973 


16.444 


0.057 


0.282 


0.895 


-0.017 


-0.094 


-0.019 


-0.019 


-0.009 




38* 


11.005338 


12.0634 


18.431 


0.309 


0.577 


0.943 


-0.004 


-0.098 


-0.003 


-0.017 


-0.011 



Notes to Table 1. 

The SMC OGLE-II identifications are given in the first column as the field number (integer part) and star number (fractional 
part). In three cases identifications in the neighboring fields were available; these are given in parentheses following the primary 
identification. The 02 — 04 data of the primary identifications have been used as criteria for inclusion in this table. Asterisks in 
the first column identify systems which have been given half weight in the calibrations of the PLC-/3 relation. The de-reddened, 
absolute-magnitude data used in the paper have been derived using: V m ax = Imax + (V — I), (V — I)o = (V — I) — Ev-i, 
My = Vmax — (m — M)o — Ay, where (m — M)o = 18.8 and Av = 2.5Ey-i, with Ey—i = 0.12. For more detail, consult the 
text. 



Table 2. OGLE-II Data for (3 Lyrae systems with P > 10 days in Baade's Window 



OGLE-I 


OGLE-II 


P(d) 


Ev-i 


h {V - I) 2 


AJi 


h 


(V-I)i 


h 


v 2 




d(kpc) 


0.036 




26.307 


0.70 






15.51 


1.82 




(17.33) 


(0.21) 


(11.8 ± 1.8) 


3.051 


45.501074 


14.243 


0.77 


16.08 1.88 


0.05 


15.98 


1.83 


16.03 


17.91 


1.11 


9.4 ± 1.2 


3.079 


46.036352 


12.313 


0.87 


16.51 2.06 


0.05 


16.43 


2.07 


16.46 


18.52 


1.49 


9.3 ± 1.2 


4.002 


46.135504 


22.672 


0.67 


14.59 1.79 


0.18 


14.33 


1.99 


14.41 


16.20 


0.44 


6.6 ±0.8 



Notes to Table 2. 

The first two columns give the star identifications consisting of the respecive OGLE field numbers (integer part) and of the number 
in the field (fractional part). The reddening estimate for d > 2kpc, Ey-i, and the maximum light values of I\ and (V — 7)i are from 
Paper III. The differences between the mean and maximum levels, A/i, have been calculated on the basis of the OGLE-I light curves. 
The values of A7i have been used to convert the average OGLE-II data to maximum light level values, J2. It has been assumed that 
the color indices (V — 1)2 are the same for the average and maximum light levels. The values of M^ 1 have been calculated using 
Eq. (3) leading to the distance estimates through: V — M^ 1 = 51ogd — 5 + 2.5£V-j. The OGLE-II data are not available for the 
system #0.036; in this case the distance estimate is based on the OGLE-I data. 



